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ABSTRACT 


Ai  model  based  on  micromechanics  for  predicting  effective  viscoelastic 
stress-strain  equations  and  microcrack  growth  in  particle-reinforced  rubber 
(or  other  relatively  soft  viscoelastic  matrix)  is  described.  Geometric 
idealization  of  the  microstructure  follows  that  of  the  composite  spheres 
assemblage  and  generalized  self-consistent  scheme  originally  used  for 
linear  elastic  composites  without  damage.  The  approach  combines  a 
perturbation  analysis  of  the  matrix,  which  becomes  mote  accurate  as  the 
particle  volume  fraction  is  increased,  with  the  Rayleigh-Ritz  energy  method 
for  predicting  mechanical  response  of  the  composite.  Results  for  linear 
elastic  behavior  with  crack  growth  are  first  obtained,  and  then  extensions 
to  linear  and  nonlinear  viscoelastic  behavior  are  discussed.  It  is  shown 
that  the  elasticity  theory  may  be  easily  extended  to  predict  mechanical 
response  of  a  viscoelastic  composite,  and  that  an  approximate  equation 
governing  microcrack  growth  is  analogous  to  one  for  an  aging  elastic 
material.  Finally,  a  limited  assessment  of  the  theory  is  made  through 
comparison  with  some  existing  effective  modulus  results  and  experimental 
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data  on  a  particle-filled  rubber. 
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1.  Introduction 


Rubber  which  is  filled  with  a  large  volume  fraction  of  hard  particles 
usually  exhibits  considerable  nonlinearity  and  strain  cate-dependent 
hysteresis  over  a  wide  range  of  temperatures  and  rates.  In  this  paper  we 
develop  a  mathematical  model  from  tnicromechanical  considerations  for 
predicting  the  growth  of  microcracks  and  their  effect  on  the  elastic  and 
viscoelastic  deformation  behavior  of  such  a  composite.  Motivation  for  the 
geometric  idealizations  and  analytical  approach  comes  from  studies  of  the 
constitutive  behavior  and  microstructure  of  composite  solid  propellant. 
However,  it  should  be  clear  that  at  least  portions  of  the  theory  may  apply  to 
many  other  composites  consisting  of  continuous  matrices  reinforced  with  a 
large  volume  fraction  of  much  stiffer  and  stronger  particles  or  fibers;  in  the 
latter  case  prediction  of  response  characteristics  for  loading  transverse  to 
the  fibers  would  be  analogous  to  that  done  here. 

Figure  1  shows  the  microstructure  of  a  typical  composite  solid  propellant 
[IJ.  In  this  case  there  is  76  volume  percent  filler  consisting  primarily  of 
ammonium  perchlorate  particles  with  some  aluininum  powder  (both  of  comparable 
stiffness)  in  a  lightly  crosslinked  rubber.  The  applied  strain  is 
approximately  30  percent,  and  is  close  to  that  producing  specimen  failure. 
Only  the  largest  particles  may  be  easily  seen  as  there  is  a  broad  distribution 
of  sizes,  primarily  in  the  10  to  lOO  micron  range.  Each  particle  is 
surrounded  by  a  strongly  attached  rubber  layer.  Typically,  beginning  at 
strains  that  are  less  than  10  percent  of  the  ultimate  value,  crack- Like  voids 
start  to  appear  near  the  largest  particles  and  grow  steadily  and  separately 
from  one  another;  a  crack  that  grows  around  one  large  particle  does  not 
usually  join  with  other  cracks  (until  the  ultimate  strain  is  approached) 
because  of  the  shielding  effect  of  adjacent  particles. 

Uniaxial  stress-strain  behavior  of  a  similar  propellant  (with  78  volume 


percent  particles)  is  shown  in  Fig.  2  over  a  two  minute  test  period.  The 
points  are  experimental  values  and  tte  continuous  line  is  from  the  equation 
developed  in  [2].  One  specimen  grip  was  allowed  to  slip  in  compression  (with 
respect  to  the  machine  crosshead)  to  avoid  buclcling  of  the  long  bar-shaped 
specimen;  the  data  points  reflect  this  low  compressive  resistance.  On  the 
basis  of  the  earlier  study  [2],  data  on  dilatation  of  propellant  [3],  and  the 
characteristic  behavior  of  stronger  rubber  with  and  without  particles  [4],  it 
is  believed  the  stress-strain  behavior  in  Fig.  2  reflects  hysteresis  due  to 
viscoelasticity  (the  difference  between  curves  2  and  3),  stiffening  due  to 
high  axial  straining  of  at  least  a  fraction  of  the  polymer  network  chains  (the 
upper  end  of  curve  3),  softening  due  to  mictocracks,  and  the  type  of  softening 
usually  called  the  Mullins'  Effect.  This  latter  phenomenon,  which  is  very 
significant  in  carbon  black-filled  rubber  without  apparent  microcracking,  is 
believed  to  be  primarily  due  to  non-affine  displacement  of  network  junctions 
and  entanglements  during  loading,  followed  by  partial  or  delayed  recovery  of 
these  points  at  reduced  strains;  with  weak  crosslinks,  more  or  less  permanent 
breakdown  of  network  junctions  may  also  occur  [4].  Figure  3  illustrates  the 
Mullins'  Effect  without  viscoelastic  hysteresis. 

In  view  of  these  several  apparently  important  mechanisms  of  inelastic 
behavior  which  influence  mechanical  response  of  the  composite,  it  is  believed 
essential  to  begin  our  study  with  a  highly  simplified  microstructural 
geometry.  When  the  effects  of  the  different  mechanisms  are  understood,  less 
idealized  geometry  may  be  introduced  as  necessary.  Accordingly,  as  a  starting 
point,  we  shall  use  a  geometry  originally  proposed  by  Kerner  (5)  for 
isotropic,  linear  elastic  media  without  cracks,  which  is  illustrated  in  Fig. 
4.  A  two-phase  composite  sphere,  consisting  of  a  spherical  particle  and 
matrix  shell,  is  embedded  in  an  effective  (homogeneous)  medium  in  Kerner's 


model.  The  bulk  and  shear  moduli  of  the  effective  medium  are  initially 
unknovm;  but  they  are  found  from  the  condition  that  the  total  composite  body 
in  Fig.  4,  consisting  of  the  two-phase  sphere  in  the  much  larger  effective 
medium,  is  to  have  the  same  moduli  as  the  effective  medium.  Hashin  [6]  calls 
this  model  the  "generalized  self-consistent  scheme".  Corrections  and 
improvements  to  this  linear  model  have  been  made  by  Smith  17]  and  Christensen 
and  Lo  [8];  additional  discussion  and  comparison  with  experimental  data  may  be 
found  in  [9].  If  the  displacements  specified  on  the  outer  boundary  of  the 
effective  continuum  are  those  for  a  uniformly  strained  homogeneous  body  in 
simple  shear,  for  example,  then  the  effective  shear  modulus  may  be  found. 
Instead,  if  displacements  for  simple  shear  ate  applied  directly  to  the  outer 
surface  of  the  matrix  shell,  then  one  obtains  the  upper  bound  to  the  shear 
modulus  for  Hashin's  composite  spheres  assemblage  [lOj.  Both  locations  for 
the  boundary  conditions  are  considered  in  this  paper,  although  expressed  in 
terms  of  applied  displacements  which  produce  axisymmetric  deformation  in  order 
to  simplify  the  analysis  when  cracks  exist. 

In  Section  2  we  propose  a  modification  of  the  geometry  in  Fig.  4  to  allow 
for  cracks  and  then  outline  the  method  of  analysis.  Sections  3-5  give  details 
of  the  approach,  which  consists  of  a  combination  of  perturbation  analysis, 
Rayleigh-Ritz  approximate  energy  method,  and  fracture  mechanics.  Rigid 
particles  and  linear  elastic,  incompressible  behavior  of  the  matrix  are 
assumed  for  most  of  the  work  in  this  paper  in  order  to  develop,  with  a  minimum 
of  complexity,  a  mathematical  model  that  has  some  of  the  important  features  of 
the  actual  composite  and  to  gain  insight  for  generalizing  the  analysis. 
Nonlinear  and  viscoelastic  effects  in  the  matrix  are  discussed  briefly  in 
Sections  6-8.  Finally,  some  numerical  results  and  comparisons  are  given  in 


The  complex  geometry  of  an  actual  particulate  composite,  such  as  that 
in  Fig.  1,  is  idealized  using  the  geometry  in  Fig.  5.  A  representative 
particle  is  assumed  to  be  sphericaL  It  is  contained  in  a  concentric  shell 
of  matrix  material,  a<  r<  b,  which  may  have  one  or  more  cracks.  This 
composite  spherical  particle  is  embedded  in  an  uncracked  shell  of  an 
effectively  homogeneous  isotropic  material,  b<  r<  c,  having  effective 
properties  of  the  particulate  composite  without  cracks.  Finally,  the  outer 
region,  r>c,  is  assumed  to  be  large  enough  to  permit  treatment  as  being 
infinite  in  extent;  it  is  considerered  to  be  effectively  homogeneous  with 
the  effective  anisotropic  properties  of  the  composite  with  particles  and 
cracks,  but  this  paper  does  not  explicitly  account  for  the  outer  region. 
The  thickness  of  the  shell  without  cracks,  b<  r  <c,  depends  on  the 
microcrack  density  in  the  total  composite,  which  is  related  to  the  size 
distribution  of  particles  Namely,  cracks  tend  to  form  and  grow  around  the 
largest  particles  in  the  real  composite  each  large  particle  with  at  least 
one  adjacent  matrix  crack  is  surrounded  locally  by  composite  material 
consisting  of  smaller  particles  in  apparently  uncracked  matrix  material;  it 
is  this  latter  material  that  is  represented  by  the  shell  b  <  r  <  c.  The 
total  composite  could  be  thought  of  as  an  assemblage  of  composite  regions 
,  r  <  c  ,  each  of  which  contains  one  large  particle  with  one  or  more 
adjacent  cracks;  if  necessary,  one  could  allow  for  a  distribution  of  sizes, 
considering  variations  in  a,b,  and  c. 

As  in  other  types  of  self-consistent  models,  the  composite  model  in 
Fig.  5  is  subjected  to  outer  boundary  displacements  at  r>>c  which 
correspond  to  those  for  a  uniform  strain  field.  The  displacement 
components  u^  (i=l,2,3)  referred  to  an  orthogonal  set  of  Cartesian 


coordinates  are,  for  the  far-field 


where  are  the  components  of  a  spacewise  constant  strain  tensor. 

(The  summation  convention  is  employed  in  which  summation  over  the  range  of 
a  repeated  index  is  implied.)  This  displacement  representation  is 
completely  general  for  small  or  large  strains,  apart  from  rigid  body 
translation  and  rotation.  The  latter  motion  may  be  included  by  adding 
constants  and  using  nonsymmetric  coefficients  in  Eq.  (1)  [11].  In 
developing  stress-strain  equations  it  is  of  course  not  necessary  to  include 
the  rigid  body  motion. 

For  simplicity  we  shall  consider  here  only  small  strain  theory  and 
displacements  which  are  symmetrical  with  respect  to  the  x-j  axis  in  Fig.  6. 
There  may  be  one  or  two  ax i symmetrical  cracks,  each  centered  at  one  of  the 
poles,  as  in  Fig.  5.  For  this  axisymmetr ical  deformation  case,  ^22  ~  ^11 
and  £^^2  “  ^13  =  ^  23  “  that  Eq.  (1)  in  terms  of  the  spherical 

coordinates  becomes 

Uj^  =  r  sin  0  cos  <t>,  ^2  =  r  sine  sin  <p,  =  r  ® 


The  displacements  in  the  radial  and  latitixiinal  directions,  respectively, 
are  found  to  be 

Uj.  =  r[ey/3+ejj(2  cos^  6  -  sin^  0)/2]  (3a) 

Uq  =  -3re^(cos  6  sin  0)/2  (3b) 


where  e^  and  e^j  are  the  dilatation  and  deviator ic  strain,  respectively, 

®v  ^  2^11  ^33  '  ®d  “  ^^^33  ■  ^11^/^ 

Use  of  the  measures  of  strain  defined  in  Eq.  (4),  rather  than  and  £33, 

is  helpful  in  the  subsequent  analysia 


In  order  to  solve  the  problem  of  Fig.  5,  one  makes  the  displacements 
in  the  outer  region  equal  those  in  Eq.  (3)  at  r>>c,  as  is  done  for  the 
problem  in  Fig.  4.  Without  cracks,  the  problems  in  Figs.  4  and  5 
obviously  are  the  same.  In  this  case,  the  exact  linear  elastic  solution 
for  the  mechanical  state  and  the  effective  Young’s  modulus  may  be  obtained 
just  as  for  simple  shear  [9]  by  a  separation  of  variables  method  in  which 
the  particle  and  each  shell  have  the  same  form  of  solution.  In  fact,  one 
finds  that  the  dependence  of  the  displacements  on  6  is  that  in  Eq.  (3) 
(which  is  from  the  uniform  strain  field)  but  functions  of  r  appear  in  place 
of  e^  and  e^j;  the  function  of  r  which  replaces  in  Eq.  (3a)  is  in  general 
different  from  that  replacing  e^  in  Ec^  (3b).  Each  shell  and  the  particle 
has  this  form  of  displacement.  An  analogous  situation  exists  for  the 
stresses. 

With  cracks  this  simplicity  certainly  does  not  exist  for  exact 
solutions.  However,  we  shall  use  it  for  developing  approximate  solutions 
based  on  the  minimum  potential  energy  principle.  For  example,  with 
reference  to  the  inner  shell  with  cracks,  a  <  r  <  b,  the  interface 
displacements  at  r=b  are  taken  as 

Uj.  =  b[q2^/3  +  q2(2  cos^  6  -  sin^  t')/2]  (5a) 

U0  =  -3b  q^Ccos  6  sin  0)/2  (5b) 

where  q^^,  q2,  q3  are  free  parameters  which  are  found  by  minimizing  the 
strain  energy  for  a  given  crack  geometry.  The  outer  material,  r  >b,  then 
deforms  as  if  there  were  an  uncracked  inner  shell  with  altered  properties. 
In  practice,  for  the  type  of  highly-filled,  rigid  particle  system  in  Fig. 
1,  the  constraining  effect  of  small,  closely  spaced  particles  just  outside 
of  the  soft  matrix  shell  around  each  large  particle  tends  to  preclude  high 
crack-induced  displacement  gradients  for  r > U  Thus,  Bq.  (5)  may  actually 


be  a  more  realistic  representation  than  would  be  the  case  if  the  material 
for  r  >  b  were  truly  homogeneous. 

The  same  form  of  displacement  variation  may  be  used  at  the  interface 
r=c,  but  with  three  additional  free  parameters,  also  to  be  found  by 
minimizing  the  strain  energy.  However,  analysis  of  the  entire  composite  in 
Fig.  5  has  not  yet  been  completed.  Rather,  a  more  limited  set  of  results 
has  been  obtained  in  which  displacements  in  Bq.  (3)  are  applied  directly  at 
r=b  for  one  case,  and  at  r=c  for  another  case.  With  these  results,  we  may 
determine  if  the  essential  characteristics  of  the  model  with  crack  growth 
agree  with  experimental  data. 

Consistent  with  the  limited  goal  of  studying  the  essential  behavior  of 
the  model,  expecially  that  due  to  the  cracked  shell,  a  <  r  <  b, 
we  assume  the  particles  are  rigid  and  the  matrix  is  incompressible.  We 
have  done  some  work,  as  yet  unpublished,  allowing  for  realistic  values  of 
rubber  compressibility  and  particle  moduli.  It  was  found 
that  only  the  effects  of  rubber  compressibility  on  the  composite  properties 
were  not  negligible  for  the  type  of  composite  in  Fig.  1;  but  they  were 
small  enough  to  be  neglected  in  a  study  of  primary  model  characteristics. 

The  effective  stress-strain  equations  of  the  composite  will  be 
obtained  by  evaluating  the  total  strain  energy  for  a  unit  volume  of  the 
composite  and  then  using  virtual  work  to  define  the  stresses, 

6W  =  a-  (6) 

where  W  is  the  strain  energy  per  unit  volume;  W  depends  on  the  microcracks, 
but  they  are  considered  to  be  of  fixed  length  in  the  virtual  deformation 
process.  For  axisymmetric  deformation  and  using  e^  and  e^  in  Eq.  (A)  as 
the  measures  of  strain,  Bq.  (6)  becomes 

5W  =  (O33  -c^3)fe^  +  (2^33  +  ‘^’33)'^ev/3  (7) 


which  yields  the  following  deviatoric  and  dilatational  equations 


^  ^33  ~  ^11  ”  »  %  =  ^^°11  ^  3W/3e^  (8) 

These  results  allow  for  material  nonlinearity.  For  an  incompressible 
composite  =  0  and  only  the  first  relation  in  Eq.  (8)  applies;  in  this 
case  Eg.  (4)  yields 

^11  "^33/2  »  =  £33  (9) 

3.  Field  Equations  and  Boundary  Conditions  for  the  Matrix  Shell 

Displacement  and  strain  distributions  in  the  matrix  shell  or  layer, 
a  <r<b,  in  Fig.  5  will  be  constructed  for  incompressible  behavior  by 
deriving  the  displacements  u^  and  u^  fron  the  curl  of  a  vector  potential 
[11].  For  syimvetry  about  the  X3  axis.  Fig.  6,  they  beccme 

Uj.  =  (r  sin  9)"^  3(ij^  sin  0)/96  (10) 

Uq  =  -(r  sin  0)"^  9(r  sin  9)/9r  (11) 

It  is  helpful  to  replace  the  potential  4;  by  a  dimensionless  displacement 
function, 

F  r  r  li-  sin  0/R^  (12) 

where  R  is  the  mean  radius  of  the  matrix  shell.  Fig.  7.  In  addition,  for 
development  of  a  perturbation  solution  we  replace  r  in  favor  of  the 
dimensionless  radial  coordinate  x,  which  is  defined  by 

r  =  (1  +  hx)R  (13) 

The  dimensionless  thickness, 

h  ^  H/R  (14) 

is  assumed  to  be  small  enough  to  use  it  as  a  perturbation  parameter. 
Observe  that  inner  and  outer  radii  of  the  matrix  shell  are  given  by  x  =  -1 
and  X  =  1,  respectively.  In  view  of  Eq.  (13), 


,  -Jb  -Ji  *,>  ’ ’ 


WJ*:*  i .  n'»t'»r»v^ 


»» 1 1  'j  »  ; 
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3/3r  =  (Rh)"^  3/3x 

(15) 

# 

V 

With  a  change  of  variables  to  F  and  x,  as  well  as  the 

use  of  the 

f 

•  ^  < 

dimensionless  radius. 

m 

»  >  , 

,v  «*• 

r'vi 

• 

A 

r  s  r/R  =  1  +  hx 

(16) 

Eqs.  (10)  and  (11)  become 

•> 

=  R  r"^  (sin  6)"^  3f/90 

(17) 

4j> 

U0=-R(h  r)"^  (sin  6)“1  3f/3x 

(18) 

f- 

The  strain-displacement  equations  [11]  in  terms  of  the  dimensionless 

•» 

A 

variables  x  and  r  are 

£j.  =  (Rh)~^  3u^/9x 

(19a) 

• 

Eq  =  (Rr)"^  (Suq/  99  +  Up) 

(19b) 

=  (Rr)"^  (Uj.  +  cote  U0) 

(19c) 

c. 

lj,0=  (Rr)"^  (9uj./39-  ug)  +  (Rh)"^3u0/9x 

(19d) 

# 

Now  rewrite  the  strains  using  the  displacements  in  Bqs.  (17)  and  (18)  , 

iC, 

Ej.  =  (h  sin  0)"^  r"2  [32F/3e3x  -  2h  r~^  3F/3e] 

(20a) 

c 

=  (h  sin  0)"^  r“2  [_32F/3e3x  +  cot  6  9F/3x  +  hr"^  3F/3e] 
u 

(20b) 

£,  =  (h  sin  0)“^  r~^  [hr“^  3F/30  -  cot  6  3F/3x) 

<P 

(20c) 

Yj.0=  (h  r)"2  (sin  0)”^  [2h  3F/3x  -  r 

V 

+  h2  r'^O^F/Se^  -  cote  3F/3e)] 

(20d) 

fc*. 

The  equilibrium  equations  [11]  for  the  radial  and 

latitudinal 

1 

directions  using  the  dimensionless  variables  are,  respectively. 

•  «  •  . 

■'I'-': 

i 

r~2  h"^  3(r^  o,.)/3x  +  (sinO)"^  3(sin0n  )/3e  -3  ■,  =  0 

L  L(; 

(21a) 

r~^  h“^  3(r^a„,.)/3x  +  (sin  fl)~^  3  (sin  (■ .-  ) /.ic*  -  cot  fi  = 

^  L  t)  ip 

^  0  (21b) 

where  o  is  the  mean  normal  stress. 


C 


a 


(22) 


E  (o^  +aQ  +a^)/3 
For  a  linear,  elastic,  isotropic  matrix,  the  stress  strain  equations 

are 

“t  2GC,  ,  0(,  =  0+  2Gej, 

where  G  is  the  shear  modulus.  The  mean  stress  o  is  not  determined  by  the 
strains  in  view  of  the  incompressibility  assumptioa  When  the  strains,  Bq. 
(20),  are  substituted  into  the  stresses,  Eq.  (23),  and  then  the  stresses 
put  into  Ec^  (21),  we  obtain  two  differential  equations  for  the  mean  stress 
a  and  displacement  function  F. 

For  a  nonlinear  elastic  material  in  which  the  strain  energy  density 
depends  on  the  shear  strain  invariant  given  in  Bq.  (61) ,  the  stress- 
strain  equations  ate  the  same  as  those  in  Eq.  (23)  except  G  is  a  function 
of  this  strain  invariant. 

Using  Eq.  (5),  the  boundary  conditions  on  the  displacements  are 
specified  at  x  =  1  as 

Uj,  =  R(l+h)  tq2^/3  +  q2  (2  cos^  6  -  sin^  G)/2]  (24a) 

U0  =  -3R(l+h)q3  (cos  9  sin  e)/2  (24b) 

Allowing  for  rigid-body  vertical  translation  of  the  particle,  q^,the 
conditions  at  x  =  -1  are 

Uj.  =  q^  cos  0  (25a) 

Ug  =  -q4  sin  0  (25b) 

4.  Perturbation  Analysis  for  a  Linear  Elastic  Matrix  Shell 

Using  a  standard  method  (12] ,  we  expand  F  and  other  solution  functions 
in  a  power  series  in  the  (assumed)  small  parameter  h, 


F  =  +  h  Fi  +  h-^Fo  +  ... 


The  relevant  functions  are  then  substituted  into  the  field  equations  and 
boundary  conditions,  and  the  equations  arising  from  each  set  of  terms  of 
the  same  order  in  h  are  solved  sequentially,  starting  with  the  lowest 
order.  in  using  this  approach,  it  is  nelpful  to  select  the  relevant 
dependent  variables  in  such  a  way  that  the  lowest  order  approximations  are 
independent  of  h.  We  find  from  an  initial  trial  of  the  technique  that  Fq 
is  indeed  independent  of  h  and  that  the  corresponding  normal  stresses  ate 
0(h"^).  (Standard  notation  is  used  here  in  that  a  quantity  of  order  in 
h  is  denoted  by  O(h^)  and,  by  definition,  0(h^)/h'^  is  bounded  as  h 0 J 


Thus,  vie  introduce  the  new  variables. 


s  =  h-^o, 


Sr  =  h^a^, 


0» 


The  modified  mean  stress  s  is  expanded  like  F  in  Bq.  (26) , 


s  =  Sq  + 


h  Sj^  +  h^S2  +  ... 


Next,  in  terms  of  the  series  for  s  and  F  express  the  equations  of 
equilibrium,  Eq.  (21),  using  Eqs.  (20),  (23),  and  (26)-(28),  and  then 
combine  the  terms  into  groups,  each  of  which  is  multiplied  by  a  common 
power  in  h.  Considering  those  with  the  h°  coefficient,  the  radial  and 
latitudinal  equations  yield,  respectively. 


9Sq/9x  =  0 


aSr/ae  -  G(sin  0)"^  a^Fg/ax^  =  O  (30) 

Surprisingly,  the  equations  coming  from  the  factor  h^  are  the  same  as  Eqs 
(29)  and  (30)  except  s^  and  replace  s  and  F^,  respectively. 


w  v.v  'aviwjvjwi  m.rjwiwi 
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The  equations  arising  from  h^  are  much  more  involved,  and  therefore  we 
shall  limit  our  analysis  to  the  zeroth  and  first  order  terms.  It  is  also 
found  that  by  neglecting  terms  0(h'^), 

®r  =  Sq  =  s^  =  s  (31) 

Inasmuch  as  the  zeroth  and  first  order  equations  of  equilibrium  have 
the  same  form,  it  is  helpful  to  work  instead  with  the  single  equation, 

ds/de  -  G(sin  6  )"^8^F/dx3  =  o  (32) 


after  noting  that  s  is  independent  of  the  dimensionless  radial  coordinate  x 
(cf.  Eq.  (29)).  This  equation  is  readily  integrated  with  respect  to  x  to 
obtain  F, 

^  -  i-  ^^2  T  ^ 

The  four  functions  of  0,  i.e.  dS/de,  f^,  fj^,  and  t2i  will  be  obtained  from 
the  four  boundary  conditions  on  displacement,  Bqs.  (24)  and  (25).  To 
accomplish  this,  it  is  helpful  to  first  express  the  conditions  on  u^  in 
terms  of  F.  Thus,  by  substituting  Eq.  (24a)  into  Eq.  (17)  and  integrating 
with  respect  to  6  we  find  for  x=l, 

F  =  (l+h)3  [-qj^(cos  0)/3  +  92 (sin  6  cos  0)/2]  +  (34) 


where  is  a  constant.  For  x  =  -1, 


,  (l-h)^  sin^  e.  . 

R  2  ^4  ii 


(35) 


and  y.2  is  another  constant.  From  Eqs.  (18)  and  (24b)  for  x=l. 


3F/3x=  3  93  h(l+h)  ^(sin^e  cos  0)/2 


(36) 
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and  from  Eqs.  (18)  and  (25b)  for  x=-l. 


9F/8x  =  h(l-h)  (sin^  9)  q^/R  (37) 

Note  that  these  conditions  on  the  latitudinal  displacement  are  0(h); 
this  implies  u  =  0  at  x  =  ±1  for  the  zeroth  order  solutioa  Substitution 
of  F,  Eq.  (33),  into  tne  four  conditions,  Eqs.  (34)-(37),  gives  equations 
for  obtaining  the  four  functions  of  9  in  Eq.  (33).  Because  the  Ug 
conditions,  Eqs.  (36)  and  (37),  are  0(h)  we  may  replace  (1+h)^  and  (1-h)  by 
unity  since  Eq.  (33)  is  only  valid  up  to  the  first  order;  consistent  with 
this  approximation,  in  Eqs.  (34)  and  (35)  we  may  use  the  approxiirations 

(l+h)3  l+3h  ,  (1-h) 2  ^  i-2h  (38) 


There  results,  finally. 


3c  /  COS  0  —  C.  n  j-  "1 

de  =  T  {  (  3-iIn  '9"  ^  ■  2  ['^2  3(q2  -q3)h  J  sin  6  cos  9 

"^4  I 

+  -J.  sin0}  (3c 


+  sin  0 
2R 


3  2  ^4  2 

f2  =  ^  *^3  h  0  0 - sin  9 


c,  -  cos  0  .  1 

=  qi(l+3h)  +  g-  [^q2+ (3q2-q3)h  Jsi 


sin  0  cos  0 


^  ^4  ,,  3,  .  2„ 

+  TC—  (h-  T)sin  0 


f  ~  cos  0  ,  r  f  l+3h .  3  1,1-2-  - 

Eq  =  -qi(l+3h)  — g —  L'^2  •  4~~ )  ~  8’'^3^J  ®  cos  9 - 


+  4  (1-Wsin^e  +  ^2  (42) 

where  the  unspecified  constants  of  integration  and  k2  have  been 
eliminated  in  favor  of  the  new  constants  c^  and  C2. 

The  shear  strain,  Eq.  (20d),  is  needed  for  the  subsequent  energy 
analysis,  and  it  is  found  to  be 
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I 


„  -2  ,  . .  r  sin  0  ds 

=  (h  r)  (sine)  [ - g—  35 


S  x+2hf^-f2l 


The  latitudinal  displacement,  Eq.  (18),  is  especially  simple  if  only  the 
zeroth  order  approximation  is  used.  For  this  case  we  find  that  it  is 
distributed  parabolically  it»  x, 

R  ds  ,,  „2. 

"  TgTT  de  ^ 

as  illustrated  in  Fig.  7. 

Also,  it  is  to  be  noted  that  the  mean  stress  a  is  derived  by 
integrating  Eq.  (39)  and  then  using  Eq.  (27),  o  =  sh"\  The  constant  c^ 
may  be  found  from  consideration  of  the  behavior  of  s;  on  the  other  hand,  C2 
has  no  effect  on  the  mechanical  variables  and  thus  may  be  taken  as  zero. 
Observe  that  there  is  a  singular  term  (sin0)~^  in  Eq.  (39),  which  cannot 
exist  in  physically  meaningful  solutions.  Removal  of  this  singularity, 
together  with  consideration  of  the  mean  stress  at  the  ends  of  the  matrix 
layer  (when  one  or  two  cracks  exist)  leads  to  the  solution  in  terms  of  q^^ 
through  q^. 

Without  cracks  we  must  set  q^^  =  0  to  remove  the  singularity  and  take 
q^  =  0  to  satisfy  the  obvious  symmetry  requirement.  Use  of  these 
conditions  in  Eq.  (39)  leads  to 

s  =  Ogh^  -  ^lq2  +  3(q2  -  q3)hlsin2  0  (45) 

where  is  an  arbitrary  constant  mean  stress.  With  only  one  crack  at  the 
bottom  of  the  particle,  cf.  Fig.  5,  set  Cj^  =  1  to  remove  the  singularity  at 
6=0  and  obtain 


•  •  ( 

■  ---.I 

t...  1 
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T  3  r  r  ciinfl  tan(0,/2)  1 

S,„ih3+  j  G  {  4  (l«h)  [in  (  |g|-  5j;^)J 

-  ^  [q2 +  3  (q2 -q2)h](sin^  e  -  sin^  0j^)  --^(cos  0  -  cos  e^)|  (46) 


where  Oj^  is  the  normal  stress  applied  to  the  matrix  surface  at  0  =  02^;  for 
an  internal  cavity  pressure  p^^,  say,  =  -p^^.  Finally,  for  two  cracks 
of  equal  size  at  the  top  and  bottom  of  the  particle,  symmetry  conditions 
imply  =  0,  and  therefore 

-  [■  [q2  +  3  (q2 -q^jh]  (sin^  6  -sin^  }  (47) 

If  the  two  cracks  are  not  of  equal  length  or  the  normal  stresses  at  the 
ends  of  a  matrix  layer  with  two  cracks  are  not  equal,  the  constant  Cj^  would 
be  used  along  with  the  constant  of  integration  in  s  to  satisfy  the  normal 
stress  condition  at  the  ends  of  the  matrix  layer.  Etesults  for  this  case  as 
well  as  for  additional  cracks  or  voids  could  be  readily  obtained,  but  they 
will  not  be  given  here. 

That  the  perturbation  solution  permits  satisfaction  of  normal  stress 
boundary  conditions  at  the  ends  of  the  matrix  layer  is  very  important  with 
incompressible  materiala  For  highly  confined  material,  such  as  the  matrix 
layer,  the  error  in  a  normal  stress  boundary  condition  does  not  decay 
within  a  few  layer  thicknesses  from  the  ends;  in  fact,  the  value  of  in 
Eqs.  (46)  and  (47)  may  have  a  strong  effect  on  the  entire  mean 
stress  distribution.  In  contrast,  a  shear  stress  condition  cannot  be 
imposed  at  the  ends  because  there  are  no  remaining  free  parameters  in  the 
perturbation  solution  for  doing  this.  Unpublished  work  of  the  author  on 
problems  for  which  an  incompressible  material  layer  is  confineo  between 
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much  stiffer  plates  shews  that  the  effect  of  a  surface  shear  stress  decays 
rapidly  and  is  relatively  small  within  a  few  layer  thicknesses  from  the 
end.  This  localized  effect  would  not  be  expected  to  significantly  effect 
the  overall  strain  energy  and  energy  release  rate  predictions  made  in  the 
next  Sectioa 

5.  Ehergy  Analysis 

The  total  strain  energy  for  the  composite  of  Fig.  5  may  be  written  in 
the  form 

+W2  (48) 

The  first  term  Oq,  6j^),  i  =  1  through  4,  is  the  strain  energy 

of  the  inner  two-phase  composite  consisting  of  a  particle  and  a  matrix 
shell,  which  can  be  evaluated  from  the  results  of  the  perturbation 
analysia  The  quantity  W2  is  the  strain  energy  of  the  shell  b<r  <c.  In  the 
present  limited  analysis  the  outer  shell  r >  c  is  omitted,  and  the 
displacements  at  the  surface  r=c  are  specified  to  be  those  of  a  strained 
continuum  using  dilatation  e^  and  deviator ic  strain  e^j  (cf.  Eq.  (3)). 
Assuming  continuity  of  displacements  Eej,  (24)  across  the  interface  r=b,  we 
may  write  W2  =  ^2^®d'  ®v'  'll'  *12'  ^  more  involved  analysis  which 

includes  wculd  not  be  needed  if  c  >>b  or  if  the  goal  were  to  obtain  an 
approximate  upper  bound  on  effective  moduli  using  the  composite  spheres 
assemblage  idealization  with  or  without  the  shell  b<  r<c. 

Inasmuch  as  tne  matrix  and  the  uncracked  outer  shell  are 


incompressible,  the  overall  dilatation  e^  and  volume  change  of  the  inner 
composite,  r<  b,  are  directly  related.  The  change  in  volume  of  the  inner 
composite  is 


m 


where  is  given  by  Eq.  (24a).  We  <±)tain 


AVj^  =  4  TTb^  qj^/3  (50) 

showing  that  q2^  is  the  dilatation  of  the  inner  composite,  as  expected. 
This  volume  change  is  also  that  of  the  complete  composite  sphere  r  <c,  so 
=  4  TTc^  e^/2  and 

=  (c/b)3  (51) 

The  free  parameters  q2,  q3  and  q^  are  selected  so  as  to  minimize  the 
total  strain  energy, 

3W  3W  aw 

fori  =  2,3,4  (52) 

These  three  conditions  and  Eq.  (51)  give,  at  least  implicitly,  the 
relationships 

^i  “  'ii^^d'  "^v'  ®o'  ^  =  1»2,3,4  (53) 

Energy  Release  Rate;  The  crack  growth,  as  defined  by  a  change  in 
0^  and  will  be  related  to  energy  release  rate,  ^  .  This  quantity  may  be 
defined  basically  as  the  mechanical  work  per  unit  of  new  crack  surface  area 
that  becomes  available  at  the  crack  tip  during  an  infinitesimal  amount  of 
growtfu  The  area  increase  used  in  this  definition  is  the  area  of  new  crack 
surface  projected  onto  the  local  crack  plane,  rather  than  all  of  the 
physical  surface  area  that  may  exist  in  the  damaged  material  around  a  crack 
tip.  For  an  elastic  material,  the  available  work  is  equal  to  decrease  in 
strain  energy  for  fixed  surface  displacements  (where  external  loads  ate 
applied.)  Considering  now  the  bottom  crack  in  Fig.  5  the  increase  in  crack 
area  is  2  sin  B  dB.  Thus, 


(54) 


=  -SW^SA  =  -(2TtR  sin  B)  3W^36 

=  (2  IT  sin  3W,p/ 36j^ 

A  similar  result  is  obviously  obtained  for  the  top  crack.  For  two  cracks 
which  are  of  equal  length  (0q  =  B)  and  which  grow  simultaneously, 

^  =  -(4  TT  r2  sin  6)“^  3V^j/3B  (55) 

where  54,p/9B  is  the  derivative  for  growth  of  both  cracks.  From  Eqs.  (48) 
(withW3  “  •  (51),  and  (52),  with  e^j  and  e^  fixed, 

3Wt  ^  9W2  \  3q .  3W^  _  3W^ 

^  (56) 

along  with  a  similar  result  for  3W.j./3eQ.  Thus,  to  evaluate  ^  it  is 
sufficient  to  consider  only  the  change  in  Wj^  for  fixed  values  of  the  four 
parameters  q^. 

The  overall  stress-strain  equations  are  those  given  in  Eq.  (8) ,  where 
it  should  be  recalled  that  the  derivatives  are  taken  with  the  crack  sizes 
fixed  and  that  W  is  the  strain  energy  per  unit  volume,  W,p/V;  here 
V  =  4itcV3.  Using  Eq.  (48)  without  W^, 

i=2 

Introduction  of  Eq.  (52)  yields 

=  v"^  3W2/3e^  (58) 

The  dilatational  equation  is  derived  similarly, 

Oy  =  V"^[(c/b)3  3W^/3qj^  +  3W2/3ey]  (59) 


3Wi 


3W2 


2)^i  ^  '^^2  1 

i/^ 


(57) 


In  differentiating  W2  it  is  supposed  that  Eq.  (51)  has  been  introduced 
first  so  that  W2  does  not  depend  on  q^. 


.  .''j 
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Use  of  the  Perturbation  Solution;  Hie  strain  energy  Wj^  is  only  for 
the  matrix  because  the  particle  is  assumed  rigid.  For  an  incompressible 
material,  the  strain  energy  density  of  the  matrix  may  be  written  in  the 


form, 


1  _  2 

wi  = 


where  y  is  the  "equivalent  shear  strain";  it  is  a  positive  quantity  that  is 
proportional  to  the  square  toot  of  the  second  invariant  of  the  strain 
tensor.  For  the  axisymmetric  problem  under  consideration, 

2 

y  =  2[{ej.-€^)^/6  +  +  (V^(|))V6  +  /4]l/2 

where  the  positive  square  root  is  to  be  used.  Observe  that  Eq.  (20) 

2  .  . 

implies  the  normal  strain  differences  squared  are  0(h  and  y  is  0(h~^). 
To  be  consistent  with  the  work  in  Section  4,  we  retain  only  the  zeroth  and 
first  order  terms  in  h  in  the  quantity  (h2y),and  obtain  the  very  simple 


result 


Y  =  Y, 


where  |*  |  denotes  absolute  value,  and 

1  2 
'^1  “  2  Y  r  0 

Thus,  only  the  shear  strain  appreciably  affects  the  strain  energy  in  a  thin 
shell.  Now,  using  Eqs.  (13),  (16),  and  (63),  as  well  as  the  approximation 
that  the  strain  energy  density  vanishes  in  the  layer  wherever  it  is 
cracked  (assuming  compressive  stresses  are  not  acting  across  the  crack 
surface)  yields, 


v'.% 


sin  0  w,  dr  d6 


=  2  J"  sin  0  r^  Wj^  dr  dO 

'"o  3 

=  ttG  R^h  f  J  sin  0  (1+hx)^  y^.^  dx  d6 


e  -1 

o 


•  ■.  ■  •I 


Substitute  Eqi  (43)  to  find 


Wf  =  j  fTG(R/h)^^  sin  6  (G  ^  ds/d0)^  d0 


Although  (ds/d0)^  contains  some  terms  O(h^)  it  is  algebraically  convenient 
to  not  explicitly  omit  them.  Using  Eq.  (39)  let  us  next  introduce  the 


quantity  y  , 


,  .  /c-  -cos  e\  , 

-G-l(ds/de)  =  2  sin'9  '  )  ^q2  +  ^  (q2^3)h]sin  0  cos  0 


Then  from  Bq.  (65) , 


01 

Wf  =  j'rTG(R/h)^  /  sin  6  d0 


which  can  be  integrated  analytically.  Also,  from  Bqs.  (54)  and  (56) , 

^  =  G  R  h-3  /3  (68) 

where  Yj^  =  Y(62^).  Recall  that  c^=l  when  6^  =  a  For  two  cracks  of  equal 
length,  9q  =  6  =  it_9^  and  =  q^  =  0;  we  also  find  from  Eqs.  (55),  (56), 
and  (67)  that  the  energy  release  rate  is  again  given  by  Eq.  (68). 

It  is  interesting  to  observe  that  Y^  is  proportional  to  the  shear 
strain  Eq.  (43),  at  the  edge  of  the  matrix  layer,  6  =  0|^,  if  the  first 
order  term  in  h,  2hfj^-f2,  is  neglected.  Consequently,  with  this 
approximation  and  Eq.  (63),  we  see  that  the  energy  release  rate,  Eq.  (68), 
is  proportional  to  the  strain  energy  density  at  0  =  6^.  This  result  is 
analogous  to  that  for  strain  energy  release  rate  in  a  long  strip  between 
rigid,  parallel  clamps,  in  which  the  strain  energy  density  is  that  in  the 
uniformly  strained  portion  of  the  strip  [13].  The  strain  energy  available 
to  drive  the  crack  in  Fig.  7  is  that  in  the  layer  just  ahead  of  the  tip. 
Our  analysis  does  not  account  for  details  on  the  scale  of  h  around  the 


t 


21 


I 


0 


!• 


!• 


crack  tip,  and  therefore  we  obtained  a  result  which  is  like  that  for  the 
long  strip  problem.  It  should  be  recalled,  however,  that  the  simple  result 
in  Eq.  (68)  is  valid  using  both  zeroth  and  first  order  terms  in  h,  and 
consequently  Y]^  is  not  limited  to  proportionality  in  the  local  shear 
strain. 

6.  Geometric  and  Material  Nonlinearities 

The  matrix  shell  analysis  made  thus  far  is  based  on  linear  theory. 
Even  if  global  or  average  strains  are  not  large,  the  local  shear  strain  is 
0(h“^) ,  Eq.  (43),  and  therefore  it  may  be  large  in  highly-filled 
composites.  However,  such  behavior  does  not  necessarily  invalidate  a 
geometrically  linear  theory  for  all  but  impractically  small  applied 
strains.  This  may  be  seen  by  first  recognizing  that  Eqs.  (17) -(20)  are 
valid  for  large  deformations  if  we  replace  displacements  and  strains  by 
velocities  and  strain  rates,  respectively,  and  consider  r  and  P  to  be 
Euler ian  coordinates  and  R  and  H  to  be  instantaneous  dimensions.  One  could 
then  der ive  displacements  by  integrating  the  velocities  with  respect  to 
time  after  rewriting  the  equations  in  terms  of  L.agrangian  coordinates  and 
the  initial  geometry.  When  the  displacements  are  small  relative  to  the 
respective  coordinates,  one  obtains  from  this  process  the  same  results  as 
now  in  Eq&  (17),  (18),  and  (20),  regardless  of  the  shear  strain  magnitude. 
These  relative  displacements  are  0(h“^)  as  are  the  normal  strains;  but  the 
shear  strain  is  0(h”^) ,  thus  permitting  us  to  make  this  distinction  between 
the  magnitudes  of  shear  strain  and  the  other  measures  of  deformation. 

The  shear-strain  dominated  mode  of  matrix  deformation  not  only  permits 
the  use  of  the  forgoing  geometrically  linear  tneory  with  large  shear 
strains,  but,  with  certain  types  of  material  nonl inearity,  also  leads  to 
mechanical  state  results  which  ate  not  much  more  involved  than  those  for  a 


fully  linear  theory.  There  is  not  space  in  this  article  to  cover  this 
topic  in  detail,  and  therefore  we  shall  consider  it  only  briefly  for  i 
power  law  nonlinear  matrix  material  characterized  by  strain  energy  density 

Wj^  =  Y^V(W-l)  (69) 

where  and  N  are  positive  constants  and  Y  is  the  equivalent  shear  strain, 
Eq.  (61);  for  N=1  Eq.  (60)  is  recovered  if  we  let  =  G.  The  stress- 
strain  equations  arising  from  Eq.  (69)  are  identical  to  those  in  Eq.  (23) 
except  G  is  not  constant  but  instead  is  given  by 

G  =  Gn  Y  (70) 

A  perturbation  analysis  similar  to  that  for  the  linear  case  may  be  enployed 
if  modified  stresses  like  those  in  Eq.  (27),  but  with  replaced  by 
are  introduced.  Neglecting  terms  0(h2)  for  these  modified  stresses,  the 
normal  stresses  are  found  as  before  to  be  equal,  Eq.  (31),  and  to  be 
independent  of  x.  Moreover,  neglecting  terms  0(h),  the  latitudinal 
equilibrium  equation  is  similar  to  Eq.  (32), 

ds/d  G  -  (sin  6)“^  3(G  9^F/3x^)/9x  =  0  (71) 

where 

G  -  Gn  (h^l  Y^gl  )N-1  (72) 

and  from  Eq.  (20d) , 

G  =  G^[(sin0  )"^i9%/9x2  ljN-1  (73) 

Equations  (71)  and  (73)  lead  to  results  which  are  analogous  to  those  for 
the  linear  case  except  N  appears  in  the  various  exponents  and  coefficients. 
When  terms  0(h)  are  not  neglected  in  the  latitudinal  equilibrian  equation, 
the  results  are  considerably  more  involved. 


It  IS  of  course  not  necessary  to  satisfy  the  equilibrium  equations 
exactly  when  one  uses  the  minimum  strain  energy  principle  to  develop 
approximate  solutions  Thus,  whether  or  not  terms  0(h)  are  negligible,  one 
could  use  the  zeroth  order  displacements  (i.e.,  those  from  Eqs.  (71)  and 
(73))  in  developing  effective  stress-strain  relations.  One  could  even  use 
the  form  of  the  displacement  distributions  derived  from  linear  theory,  but 
we  have  found  that  it  is  not  satisfactory  when  strong  material  nonlinearity 
exists. 

7.  Viscoelastic  Behavior 

All  of  the  results  obtained  so  far  in  this  paper  may  be  readily 
extended  to  linear  and  a  certain  type  of  nonlinear  viscoelastic  behavior. 
The  linear  viscoelasticity  counterpart  of  the  radial  stress  in  Eq.  (23), 
for  example,  is 


Oj.  =  a  +  2 


I 


G(t-t',t)  gp-dt’ 


(74) 


o 

where  G(t-t',t)  is  the  so-called  shear  relaxation  modulus,  which  imparts 
hereditary  characteristics  to  the  deformation  behavior.  The  second 
argument  in  G  (i.e.  t)  allows  for  aging;  this  aging  may  be  due  to  chemical 
or  physical  causes,  including  transient  temperature.  It  is  assumed  the 
body  is  undeformed  for  t  <0,  but  the  lower  limit  in  Eq.  (74)  and  succeeding 
hereditary  integrals  should  be  interpreted  as  0~  to  allow  for  the 
possibility  of  a  discontinuous  change  in  strain  or  displacement  at  t=Q. 

Viscoelastic  solutions  will  be  derived  from  elastic  solutions  through 
a  correspondence  principle  designated  as  CP-II  in  [14].  This  principle 
uses  time-dependent  solutions  rather  than  Laplace  transforms,  and  is 


applicable  to  compressible  or  incompressible  bodies  with  stationary  and 
growing  cracks  and  with  large  deformations.  With  large  deformations,  the 


coordinates  r  and  9  and  all  dimensions,  such  as  R  and  h,  should  be 
interpreted  as  quantities  which  refer  to  the  undeformed  geometry.  In  order 
to  indicate  how  this  correspondence  principle  would  be  used  here,  it  is 
helpful  to  introduce  a  so-called  pseudo  variable. 


z. 

f^  -  G~^J  G(t-t',t)  ||rdt' 


(75) 


where  f  represents  a  displacement- like  quantity  such  as  u_  and  u^,  or  strain 
or  generalized  displacement  q^^.  with  the  R  superscript,  a  quantity  is 
called  a  pseudo  displacement  or  pseudo  strain.  The  coefficient  Gp,  the 
"reference  modulus",  is  an  arbitrary  constant;  it  is  usually  selected  to 


have  the  dimensions  of  modulus  in  order  for  f^  and  f  to  have  the  same 


dimensiona  The  inverse  of  Eq.  (75)  is 
t 

f  = 


■%/ 


3f 


J(t-t',t)|^  df 


(76) 


where  J  is  the  shear  creep  compliance;  it  is  related  to  G  through  a 
hereditary  integral  [14].  With  this  notation,  Eq.  (74)  becomes  an  elastic¬ 
like  equation. 


Oj.  =  a  +  2Gr  e  ^ 


(77) 


A  type  of  nonlinear  viscoelastic  behavior  is  characterized  through  the  use 

of  the  strain  energy  density  in  Eq.  (69) ,  but  with  y  expressed  in  terms  of 
R  R 

pseudo  strains  ,  Y^q,  etc.,  instead  of  the  original  physical  variablea 
For  further  discussion  of  this  nonlinear  characterization  (which  includes 
nonlinear  viscous  and  elastic  behavior  as  special  cases)  see  [14-16]. 

Let  us  now  explicitly  extend  the  micromechanics  analysis  to 
viscoelasticity,  beginning  with  Eq.  (1),  in  which  the  strains  are 
considered  to  be  specified  functions  of  time.  Operate  on  Eq.  (1)  with  the 
integral  in  Eq.  (75)  to  obtain, 


According  to  CP-II,  all  of  the  earlier  analysis  is  valid  for  a  viscoelastic 
material  if  ail  displacements  and  strains  are  interpreted  as  pseudo 
variables;  i.e. ,  introduce  the  superscript  R  and  use  Gpj  for  shear  modulus 
(instead  of  G)  as  illustrated  in  Eq.  (77).  The  viscoelastic  stresses 
predicted  by  this  method  do  not  have  to  be  calculated  fron  pseudo 
stresses;  instead,  they  are  sirtply  those  in  the  elastic-like  analysis.  If 
boundary  tractions  instead  of  the  are  specified,  then  one  would  solve 
for  from  the  elastic-like  equations  in  terms  of  the  given  tractions, 
after  which  the  time-dependent  strains  may  be  calculated  as  indicated  in 
Eq.  (76). 

8.  Prediction  of  Crack  Growth 

Besides  giving  a  procedure  for  obtaining  viscoelastic  solutions  from 

elastic  solutions,  the  theory  in  [141  provides  a  method  for  calculating 

initiation  and  continuation  of  crack  growth.  We  start  here  with  the  energy 

release  rate  analysis  in  Section  5;  Eq.  (68)  will  be  used,  in  which  all 

displacement-like  quantities  qj^  in  y  ,  Eq.  (66),  are  to  be  considered  as 

R 

pseudo  variables,  qj^  .  For  simplicity,  it  is  assumed  either  one  crack  or 
two  equal-length  cracks  are  in  the  matrix  layer. 

It  may  be  helpful  to  discuss  first  crack  growth  in  a  linear  elastic 
material  in  which  crack  growth  occurs  when  a  constant  "critical"  value  of 
energy  release  rate, is  reached.  The  equation  for  predicting  the 
instantaneous  length  of  one  or  two  equal  cracks,  as  defined  by  6  = -n  -6^, 
Fig.  5,  follows  directly  from  Eq.  (68)  after  setting 

=  G  R  h"^  (79) 

A 

where,  as  before,  y  is  defined  by  Eq.  (66)  in  which  0=0^.  in  studies 
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conducted  to  date,  such  as  those  discussed  in  Section  9,  we  have  found  that 
the  crack  growth  is  stable.  That  is,  as  the  applied  strain  is  increased 
(which  causes  changes  in  qj^)  Eq.  (79)  yields  a  unique,  increasing  crack 
angle  3.  This  crack  growth,  in  turn,  yields  decreasing  effective  shear  and 
bulk  moduli,  as  determined  from  Bqs.  (58)  and  (59). 

The  crack  growth  theory  in  [14  j  expresses  the  growth  in  terms  of  a 
generalized  J  integral,  designated  as  J^,  For  so-called  "self-similar" 
growth  [14]  (which  is  concerned  with  deformation  details  on  the  scale  of 
the  crack  tip)  is  equal  to  ^  if  pseudo-variables,  rather  than  the 
physical  variables,  are  used.  Within  the  context  of  the  perturbation 
analysis,  the  value  of  ^  is  not  sensitive  to  the  difference  between  self¬ 
similar  and  more  general  growth  and  consequently  we  shall  assume  J^=*S. 
Experience  with  the  type  of  composite  in  Fig.  1  indicates  that  the 
instantaneous  crack  speed  obeys  a  power  law  in  For  a  linear 

viscoelastic  material  with  a  single  mode  of  crack  tip  deformation,  ^  is 
proportional  to  the  square  of  stress  intensity  factor;  the  power  law 
relationship  between  speed  of  maccocrack  growth  and  this  factor  is 
illustrated  in  [17]  for  solid  propellant.  Thus,  it  is  supposed  here  that 

de/dt  =  (^/^Q)Vtg  (80) 

where  is  a  positive  constant  with  dimensions  of  energy  release  rate,  and 
k  is  a  dimensionless,  positive  constant;  the  parameter  tg  is  a  positive 
quantity  (with  dimension  of  time)  which  may  vary  in  time  to  allow  for 
aging,  including  transient  temperature  effects.  One  may  also  introduce  a 
crack  growth  initiation  condition,  but  experience  with  solid  propellant 
shews  that  the  initiation  time  is  often  negligible  corpared  to  the  total 
period  of  crack  growth. 

As  an  aid  to  examining  the  effect  of  applied  strain  history  on  crack 


growth,  it  is  helpful  to  consider  here  only  those  situations  for  which 
there  is  one  independently  applied  strain  e,  and  e  and  the  associated 
pseudo  strain  e  are  non-negative.  For  example,  specify  a  time-dependent 
axial  strain  e  =  —0  ^nd  assume  the  other  components,  and  ^22'  vary 

proportionally,  or  are  such  that  a  state  of  average  uniaxial  tensile  stress 
exists  in  the  x^  direction.  Although  other  less  restrictive  cases  have 
been  studied,  we  shall  limit  our  analysis  to  the  aforementioned  situations. 
For  each  crack  length  the  shear  strain  parameter  y  ,  Eq.  (66),  is 
proportional  to  e  because  the  problem  is  linear.  Thus,  we  may  write  for  an 
elastic  material, 

=  cf(B)  (81) 

The  function  f  =  f(6)  is  derived  by  minimizing  the  strain  energy  with 
respect  to  the  free  parameters,  as  discussed  in  Section  5,  and  by  using  the 
specified  transverse  straining  or  loading  conditions,  such  as  ^ii=^22~^' 
In  all  cases  studied  to  date  we  have  found  that  f  is  a  positive, 
continuous,  decreasing  function  of  P.  Thus,  for  each  c  the  elastic  energy 
release  rate  from  Eqs.  (68)  and  (81) , 

^  =  G  R  h"3f2  e2/3  (82) 

is  a  continuous,  decreasing  function  of  6.  if  we  assume  crack  growth 
occurs  when  ^  (assumed  constant)  and  solve  Eq.  (82)  for  B  as  a 

function  of  e ,  we  obviously  would  find  that  the  microcracks  grow  stably  as 
the  applied  strain  is  increased. 

For  a  viscoelastic  material,  as  discussed  in  Section  7,  one  replaces  G 

by  Gp  and  e  by  in  Eq.  (82) ,  where 

t 

R  -1  /*  de 

^  =  Gr  /  G(t-t’,t  )dt^clf  (83) 

o 

Assuming  crack  growth  is  defined  by  Bq.  (80) , 

dp/dt  =  (e^)P  g/tg 


(84) 


where  p  ::  2k  and 


g  =  g(6)  =  (Gr  R 

Let  us  next  introduce  a  "reduced  time". 


(85) 


t 

t  =  t(t)  =  j  dtVtg(f) 
o 

A 

Without  aging  tg  is  constant  and  the  reduced  time  t 
Integration  of  Eq.  (84)  yields 


(86) 

becomes  t/t^. 


6 

J  g“^  d8  = 

where  S  is  the  Lebesgue  norm  of  pseudo-strain. 


(87) 


S  =  S(t)  E 


(88) 


o 

and  6q  is  the  initial  value  of  6.  Also,  =  c^(t')  in  Eq.  (88),  where  t' 
is  the  dummy  variable  of  integration. 

Equation  (87)  may  be  solved  implicitly  to  obtain  the  instantaneous 
crack  angle  6  as  a  function  of  the  Lebesgue  norm  S  and  6^.  inasmuch  as  the 
value  of  S  determines  the  amount  of  crack  growth,  it  will  be  termed  a 
"damage  parameter".  Because  p  is  often  quite  large,  it  is  found  that 
numerical  and  approximate  analytical  predictions  of  B  are  aided  by  using 
the  pth  root  of  the  integral  in  Eq.  (88)  for  the  damage  parameter,  rather 
than  SP  as  in  Eq.  (87).  It  should  be  noted  that  when  p=“'  ,  S  becomes  the 
largest  value  of  over  all  time  up  to  the  present  [18J.  For  values  pji  4 
and  de^/dt^.  0  an  excellent  approximation  is  12), 

S  =  eR  £l/p/(p  +  i^)l/P  (89) 

where  m  =  dlog  e  /dlog  t.  Similar  but  sonewhat  more  involved  approxima- 

p 

tions  may  be  constructed  for  histories  e  which  cure  not  monotone  increasing 
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In  many  cases  the  denominator  in  Eq.  (89)  is  essentially  constant 
because  p  is  large  or  m  does  not  vary  significantly.  If  so,  the  damage 
parameter  depends  on  only  the  current  value  of  pseudo  strain  and  time. 
This  behavior  enables  us  to  rewrite  the  basic  crack  growth  equation  so  that 
it  looks  like  one  for  an  aging  elastic  material.  Specifically,  substitute 
Eqs.  (85)  and  (89)  into  Eq.  (87),  and  rewrite  it  to  obtain. 


R  h-3(Y,^)  /3 


(90) 


where 


-  R 


f  e 


(91) 


r  ^  T  i/k 

^CR  -  %  |_  f  (92) 

and  c  is  the  denominator  in  Eq.  (89);  recall  that  p=2k.  The  tight  side  of 
Eq.  (90)  has  the  same  form  as  Eq.  (79)  for  energy  release  rate  in  an 
elastic  material.  The  quanti ty tnus  plays  the  role  of  a  critical 
energy  release  rate;  however,  instead  of  being  constant,  it  depends  on 
crack  length  through  8  and  6^  and  on  time.  It  was  observed  in  Section  7 
that  the  viscoelastic  constitutive  equations  in  terms  of  pseudo  strains 
are  analogous  to  those  for  an  elastic  material.  Equation  (90)  shows  that 
this  observation  extends  to  crack  growth  behavior  as  long  as  one  accounts 
for  crack  size  and  possibly  aging  in  the  critical  energy  release  rate. 

This  elastic-like  behavior  can  be  shown  to  exist  when  there  is  more 
than  one  independently  specified  strain  or  stresa  Indeed,  using  the  self- 
consistent  model  described  in  this  paper,  it  leads  to  the  existence  of  a 
strain  energy-like  potential  with  crack  growth.  Such  behavior  is  discussed 
in  [19],  especially  for  elastic  materials  with  damage,  and  an  illustration 
is  given  using  experimental  data  on  solid  propellant  under  confining  pressure 
and  axial  stretching. 


c 


Vf  JLV-  .r,  . 


t 


Finally,  it  should  be  remarked  that  there  ate  realistic  situations  in 
which  the  crack  growth  equation  is  that  in  F,q.  (80),  but  the  constitutive 
behavior  is  that  for  an  elastic  material  because  the  shear  relaxation  modulus 
is  essentially  constant.  Alternatively,  the  crack  growth  equation  may  be 
like  that  for  an  elastic  material,  constant,  but  the  relaxation  modulus 
is  not  constant;  this  case  may  be  recovenad  from  Eq.  (92)  by  letting  k 
9.  Some  Results  and  Comparisons 

Range  of  Validity  of  the  Perturbation  Solution;  i-i  i tnout  cracks  the 
accuracy  of  the  perturbation  analysis  in  predicting  the  effective  Young's 
modulus  of  the  composite  may  be  assessed  through  comparison  with  certain 
exact  results.  We  shall  consider  two  cases;  for  both,  it  is  assumed  the 
particle  volume  fraction  for  the  composite  sphere  of  radius  b  in  Fig.  4  is 
the  same  as  that  for  the  total  composite.  In  the  first  case  the 
displacements  in  Eq.  (3)  are  applied  directly  on  the  surface  at  r=b  in  Fig. 

4,  giving  a  modulus,  denoted  by  which  may  be  interpreted  as  an 

approximate  upper  bound  to  the  modulus  of  Hashi n's  composite  spheres 
asssnblage  110).  The  exact  expression  for  this  upper  bound  will  be  denoted 
by  E^g;  it  is  three  times  the  shear  modulus  derived  in  [10)  in  view  of  the 
incompressibility  condition. 

TO  derive  Ej^  we  first  note  that  Eq.  (24)  must  be  the  same  as  Eq.  (3); 
thus  q]^=ey  and  q2“‘l3“^d-  Moreover,  without  cracks  e.^=q4=0,  and  ejj=E23  from 
Eq.  (9).  Equation  (66)  reduces  to 

9=3  C33  sin26/8  (93) 

Substitute  Eq.  (93)  into  (67),  with  0 ^=0  and  0  ^=n,  divi.de  the  strain  energy 
W3  by  the  volume  4‘  b^/3,  recall  th.iL  R  is  the  mi-.m  i.kIius  (cf.  Fig.  /)  ,  and 


then  use  this  result  for  W  in  the  first  relation  in  Eq.  (3)  with  033-O; 
there  results,  finally. 


Ej3  =  EV'[20n^(l+3h)  J 


(95) 


is  the  composite  effective  Young's  modulus  in  terms  of  the  matrix  Young's 

modulus,  E=3G.  The  third  column  in  Table  I  gives  the  ratio  E./E  for  a 

^  b'  cs 

range  of  particle  volume  fractions  Vp;  the  second  column  comes  from  the 
exact  relationship 

h  =  (1  -  Vp  ^/3)/(l  +  Vp  1/^)  (96) 

which  may  be  easily  derived  using  the  definition  for  h  in  Eq.  (14).  It  is 
encouraging  that  Ej^  agrees  so  well  with  the  exact  result  E^g  for  particle 
volume  fractions  as  low  as  0.4. 

The  second  comparison  we  shall  make  is  for  the  generalized  self- 
consistent  model  in  Fig.  4,  in  which  boundary  displacements  in  Bq.  (3)  ate 
applied  at  a  radius  r=c  >>b.  In  this  case  the  exact  solution  for  Young's 
modulus, denoted  by  Eg^,,  may  be  obtained  from  three  times  the  shear  modulus 
in  [9,  pi  56].  The  approximate  malulus,  say  E^,  follows  from  Bq.  (53)  with 
"d  ~  °33'  ‘^d  "  ^33'  ^4  “  ~  ^  strain  energy  of  the  effective 

medium,  t>  b,  can  be  evaluated  using  exact  displacement  distributions 
similar  to  hose  in  [9].  Tlie  process  is  straightforward  and  the  amount  of 
algebra  is  considerable.  Therefore,  only  numerical  values  for  the  ratio 
Ec  given  here;  they  are  listen  in  the  fourth  column  of  Table  1. 

The  error  in  the  modulus  which  uses  the  pierturbation  solution  is  within  13'fe 
of  the  exact  solution  for  1  •  Vp  •  0.65;  for  the  composite  in  [2]  the  error 
IS  practically  zero  since  Vp  =  a78.  The  last  column  in  Table  1  gives  the 
ratio  of  exact  values,  showing  that  the  effective  modulus  of  the 
generalized  self-consistent  scheme  is  considerably  smaller  than  Hashin's 
upper  bound  for  tt>e  Vp  range  indicated;  however,  it  is  found  tnat  ‘  E^g 
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dS  Vp->  0. 

The  Perturbation  Solution  with  One  Crack;  Consider  again  the  case  in 
which  the  displacements  in  Eq.  (3)  are  applied  on  the  surface  r=b,  where 
^23^0  and  =  o  ace  the  given  inputs.  We  now  use  both  expressions  in 
Eq.  (8)  to  derive  ®v*  strain  energy  v'J  is  again  that  in  Eq. 

(67),  after  dividing  by  the  volume.  For  Y  ,  Eq.  (66),  we  use  q^=ey, 
*32"^3“^d'  ^i"l»  minimize  the  strain  energy  with  respect  to  to  find 
this  particle  displacement.  The  crack  angle  B  is  derived  from  Eq.  (79) 
assuming  is  constant. 

Figures  8  and  9  show  representative  distributions  of  the  shear  strain 
measure  Y  I  Eq.  (66),  and  the  dimensionless  normal  stress  s/G  for  three 
values  of  the  axial  strain  c  i  ^33?  recall  that  s  =  h^o  ,  Eq.  (27),  where 

A 

0  is  the  normal  stress  in  the  matrix.  In  Fig.  8  the  peak  values  of  >  are 
the  same  since is  constant.  These  plots  are  for  the  given  values  of 
Y3  =  0.09  and  Vp  =  0.78  (but  the  results  ate  not  very  sensitive  to  Vp) ;  for 
the  three  strains  the  predicted  values  of  B  ace  approximately  4,  13,  and  23 


degrees. 

Results  without  a  crack  ace  also  drawn  in  Figs.  8  and  9.  They  are 
from  Eq.  (46)  (in  which  qj^=q4=0,  q2=q3=  ,  and  ")  and  from  Eq.  (93). 

Without  a  crack  the  normal  stress  at  the  poles,  in  Lq.  (46),  is  not 
given  by  a  boundary  condition;  since  the  material  is  incompressible  the 
strains  determine  the  normal  stress  only  to  within  a  constant. 
This  constant  may  be  found  by  requiring  the  total  radial  force/sucface  area 
acting  on  the  matrix  shell  to  be  the  mean  applied  stress,  '  33'3.  v-Jith  this 
condition  we  obtain 

s/G  =  [1  +  0.  2./(l+3h)  -  1.5  sin2  Olr/4  (97) 


rl 
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Additional  numerical  results  confirm  that  the  quantity  s/G  from  Eq.  (46) 
for  one  crack  and  Oj^=o  approaches  Eq.  (97)  with  decreasing  crack  length, 
except  near  the  crack  tip.  Figures  8  and  9  show  that  the  effect  of  the 
crack  on  y  and  s/G  for  c  =  Cl05  is  not  great  except  neat  the  tip,  and  that 
for  all  three  strains  the  only  qualitative  difference  in  the  distributions 
for  no  crack,  Eqs.  (93)  and  (97),  and  for  one  crack  is  neat  the  crack  tip. 
In  additional  unpublishel  studies,  the  maximum  6  of  23  degrees  in  these 
figures  (for  c  =  0.15)  was  found  to  be  somewhat  larger  than  needed  to 
achieve  agreement  with  experimental  results  in  (21  on  stress-strain 
behavior;  thus,  in  this  case,  curves  for  the  largest  strain  exaggerate  the 
change  in  the  distributions  due  to  crack  growth. 

Studies  with  two  cracks  have  been  made  as  well.  However,  lor  the  same 
values  of  and  initial  crack  length  as  with  one  crack  it  is  found  that  a 
larger  value  of  applied  strain  is  needed  to  initiate  simultaneous 
propagation.  It  is  thus  likely  that  even  if  two  diametrically  opposed 
cracks  were  to  exist  initiaily,growth  of  only  one  would  be  significant. 
Comparisoqs  of  predicted  and  measured  stress-strain  behavior  made  to-date 
also  support  the  one-crack  model.  Predicteti  stress-strain  curves  and 
various  comparisons  will  be  published  elsewhere. 

Effects  of  Viscoelasticity;  Incorporation  of  viscoelasticity  into  the 
model  has  been  discussed  in  Section  7.  Essential  aspects  of  the  method  for 
dealing  with  both  crack-tip  and  global  effects  are  confirmed  by  the 
experimental  results  in  [2].  Here  we  shall  discuss  only  the  theoretically 
derived  damage  parameter,  Eq.  (88),  and  the  relationship  of  this  Lebesgue 
norm  to  the  findings  in  [2J.  Neglecting  aging  effects,  the  reduced  time  t, 

A 

Eq.  (86),  reduces  to  t  =  t/tg  in  which  tg  is  a  constant.  The  exponent  p  is 
related  to  the  exponent  which  character izes  time-dep»jna..}nce  of  the  stress 


relaxation  modulus  and  to  the  characteristics  of  the  material  in  the  zon? 
of  failing  material  at  a  crack  tip.  For  a  nonaging  material  G(t-t',t,' 
reduces  to  G(t-t'),  and  thus  for  power  law  relaxation, 

G(t)  =G;,^t"‘^  (98) 

where  Gj^  and  n  are  positive  constants.  According  to  the  data  for  rubber  in 
[20: III],  the  failure  zone  at  the  crack  tip  is  elastic-like,  which  in  turn 
yields  [20:Il], 

p  =  2(1  +  1/n)  (99) 

This  relationship  is  in  agreement  with  data  for  mactocrack  growth  in 
solid  propellant  [17].  For  the  relaxation  modulus  in  [2]  n  =  0.366,  and 
thus  p  =  7.46.  From  the  same  study  the  exponent  was  found  experimentally 
to  be  p  =  6-5.  The  difference  may  be  in  part  due  to  the  material 
nonlinearity;  namely,  apart  from  effects  of  damage,  a  power  law 
nonlinearity  was  reported,  in  which  'v  (t:^)*^,  where  N  =  0.85.  In 
another  micromechanics  model  study  it  was  shown  that  instead  of  Eq.  (99) 
one  should  use  p  =  (1+N)k  [16,  Eq.  (142));  the  quantity  k  is  the  same  as 
for  a  linear  viscoelastic  material  [15],  k  =  (1+1/n).  Thus,  p  =  6.9,  which 
is  only  6%  greater  than  the  experimental  value  of  6.5. 

The  experimental  data  in  [2]  indicate  that  two  different  Lebesgue 
norms  are  needed,  corresponding  to  p  =  6.5  and  p  =  ,  where  the  latter 

norm  is  equal  to  the  largest  value  of  T^^lt')  over  all  O^t'^t,  where  t  is 
the  current  time  [18].  The  value  p  =  «>  is  what  one  would  predict  for  an 
elastic  material  since  n  =  0.  The  power  law  in  Eq.  (93)  with  n  =  0.366  is 
for  a  broad,  intermediate  time  range,  where  in  crack  growth  theory  t  is 
one-third  the  time  taken  by  a  crack  in  growing  an  amount  equal  to  tn..- 
length  of  its  crack  tip  failure  zone  [20:11].  Thus,  considering  that  there 
is  actually  a  distribution  of  particle  and  microcrack  sizes,  and  associate! 


microcrack  growth  rates,  the  need  for  p  =  °"  may  arise  from  tiie  very  slow 
and  fast  growth  rates  for  which  elastic-like  behavior  would  exist. 
Additionally,  molecular  scale  damage  processes  related  to  breakdown  of 
crosslinks  and  highly  entangled  long  molecular  chains  (normally  associated 
with  the  so-called  Mullins'  effect  [41)  may  account,  at  least  in  part,  for 
the  need  of  the  p  =  “  value.  The  value  of  p  may  also  depend  on  whether  or 
not  the  crack  tip  is  in  the  matrix  or  at  the  particle-matrix  interface;  in 
this  regard  it  should  be  noted  that  the  perturbation  results  are  not 
sensitive  to  the  radial  location  of  the  cracks,  and  thus  they  apply  to  both 
cohesive  and  adhesive  fracture  processes.  Although  it  is  not  yet  clear 
just  why  p  =  =”  is  needed,  it  is  very  encouraging  that  the  two  values  of 
p  =  6.5  and  p  =  »  bracket  the  theoretically  predicted  range  of  possible 
exponents  using  viscoelastic  fracture  mechanics  for  the  matrix. 

With  exponents  of  p  =  6.5  and  larger,  cottespondinj  to  k'!^  3.2,  we  have 
found  that  the  value  of  the  effective  critical  fracture  energy,  Eq.  (92), 
changes  only  a  small  amount  with  straining  except  for  t  very  close  to  its 
initial  value  Indeed,  curves  analogous  to  those  in  Figs.  8  and  9,  but 

for  the  viscoelastic  material,  are  practically  the  same  as  for  elasti: 


behavior. 


It  is  also  of  interest  to  observe  that  for  t  =  t/t^  and 


o  rn 

e  =Cj^(t/t  )  ,  with  Cj^  and  m  constant  (such  as  would  be  obtained  for  a  pwer 


law  relaxation  modulus  and  a  constant  strain  rate  input)  we  may  eliminate 
frcttv  Fjq.  (92)  using 


/S 

t  =  (iOO) 

When  this  is  done  the  resulting  expressions  are  the  same  as  in  Eqs.  (93) 
and  (92)  except  the  exponent  on  is  2  +  1/mk,  is  .aultiplied  by 
,  and  '^cR  no  longer  depends  explicitly  on  t.  Namely,  the  crack 
growth  behavior  is  pr^xlicted  to  be  the  same  as  for  a  nonl  inear  powei  1  iw 


n  > 


y 

I 

L'  -• 
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elastic  material  with  a  strain  rate  and  crack  Lenytti- iep-^ndent  fractui 
energy  =  ^cR^^l  ,6).  Also,  as  a  result  of  the  incegtal  in  Eq.  (92;, 
this  energy  is  initially  zero  and  then  initially  itu  census  wit:i  .•tac-if 
growth. 

10.  Conclusions 

A  model  has  been  described  for  predicting  mictoccnck  growth  ami  tne 
overall  or  effective  constitutive  equations  of  a  nighly-f i 1 Imi  particuljte 
composite.  Assuming  rigid  spherical  particles  an^i  an  incompcessiole 
matrix,  the  theory  was  developed  using  a  combination  of  jiectucbat ion  ami 
strain  energy  methods.  Accuracy  of  the  perturbation  analysis,  whicn 
increases  with  increasing  values  of  particle  volume  fraction,  Vp,  was 
assessed  by  comparing  the  effective  Young's  modulus  with  existing  exact 
values  predicted  for  the  composite  spheres  assemblage  and  the  generalized 
self-consistent  model.  In  the  latter  case  the  error  is  less  than  13  %  for 
Vp>0.65,  and  in  the  former  case  the  error  is  less  than  8%  for  Vp'-0.40. 
Viscoelastic  effects  were  introduced,  arxl  a  damage  parameter  that  defines 
the  extent  of  time-dependent  microciracking  was  derived  and  found  to  Cue  in 
good  agreement  with  experimental  results  on  a  rubber  filled  with  relatively 
rigid  particles. 

This  paper  described  only  the  basic  portion  of  a  much  more  extensive 
study  underway  to  predict  the  mechanical  response  of  linear  and  nonlinear 
viscoelastic  particulate  composites.  The  geometric  idealization  of  the 
composite  and  the  method  of  analysis  were  selected  to  be  realistic  but  yet 
simple  enough  that  it  would  be  feasible  to  includ'i  geometric  and  material 
nonl  inear  ities,  compressibility  of  the  matrix,  and  parti'^le  deformation 
under  multiaxial  loading  of  the  composite.  Studies  beyond  those  coveretl 
here  have  been  made  and  will  be  reported  elsewhere. 
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uniaxial  stress  (in  the  direction  of  the  arrow)  just 
below  that  causing  specimen  failure,  as  observed  through 
a  scanning  electron  microscope  (lOOX).  After  Cornwell 
■  and  Schapery  [l], 
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Illustraclon  of  the  >ftilllns*  Effect  in  filled  rubber  without 
viscoelasticity  and  microcracking. 


Figure  4.  Generalized  self-consistent  scheme  without  raicrocracking 
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